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PREFACE 


The progress on the research carried out under NASA Grant NsG 3012 is 
reported briefly here. The period covered is for six months ending December 31, 
1977. 


B . Lakshrainarayana 
Principal Investigator 
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1. STATEMENT OF THE PROBLEM 


Both the fluid dynamic properties of a compressor rotor, governing efficiency, 
and the acoustical properties, governing the radiated noise, are dependent upon 
the character of the wake behind the rotor. There is a mixing of the wake and 
the free stream and a consequential dissipation of energy, known as mixing loss. 

The operational rotors generate the radial component of velocity as well as the 
tangential and axial components and operate in a centrifugal force field. Hence, 
a three-dimensional treatment is needed. 

The objective of this research is to develop an analytical model for the 
expressed purpose of learning how rotor flow and blade parameters and turbulence 
properties such as energy, velocity correlations, and length scale affect the 
rotor wake charadteristics and its diffusion properties. The model will necessarily 
include three-dimensional attributes. The approach is to employ, as information 
for the model, experimental measurements and instantaneous velocities; and 
turbulence properties at various stations downstream from a rotor. A triaxial 
probe and a rotating conventional probe, which is mounted on a traverse gear 
operated by two step motors, will be used for these measurements. The experimental 
program would include the measurement of mean velocities, turbulence quantities 
across the wake at various radial locations and downstream stations. The ultimate 
objective is to provide a rotor wake model, based on theoretical analysis and 
experimental measurements, which the acousticians could use in predicting the 
discrete as well as broadband noise generated in a fan rotor. This investigation 
will be useful to turbomachinery aerodynamists in evaluating the aerodynamic 
losses, efficiency and optimum spacing between a rotor and stator in turbomachinery. 
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2. THEORETICAL INITISTIGATION 

The material presented here is a continuation of the analysis described 
in reference* [1] . 


2.1 Analysis Based on k-g Model 

2.1.1 Modification of turblence dissipation equation to include the effects of 
rotation and streamline curvature 

To describe the turbulence decay, second order dissipation tensor** = 


u(3u^/3x^) (3Uj/3Xj^) should be modeled. For simplicity of calculation and modeling, 
the scalar representation of dissipation is usually employed in representing the 
turbulent flow field. Scalar dissipation equation can be derived through the 
contraction of dissipation tensor e = stationary generalized coordinate 

system, the transport equation for e can be written as follows for high turbulent 
Re 3 molds number flows by retaining only higher order terms [2,3] 


I 


r 1 , ,ik j . , 3,k 1 , , j,i k . 

-v[s* u,,u. . + s’ U, ,u, . + s!,u-” U, . + si, u-” u,.] 

k 1,J 1 k,j ik ’j ik j 

II 


+ P 1 - Ceu^) 


III 


IV 


2 mlt r T " ■ t T 1 


( 1 ) 


Term I of the above equation is lower order for flows with straight streamline 
but is larger for flows with streamline curvature and hence it is retained. 

Jn a rotating frame, additional terms due to centrifugal and Coriolis forces 
enter into the momentum, dissipation equations. The additional term which appear 
in the right hand side of Eq. 1 is as follows [5] 


*Numbers in square brackets denote the references listed at the end of this report. 

**A list of notations is given in Appenidx II. 
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+ £, s‘ u, . 
ipq ’k kpq i 


iik q ] 


’• Jlpq ’m mpq ’JJ-" 

3 

The component in Cartesian frame turns out as follows 


,„3, 3u 2 SUj Su^ 8^3 a„3 SU2 

-^Q, [- — — — + -r— — + 


3x^ 3X2 


9x2 ' 9x^ 9X2 


9uo 3u, 

1 i] 

9X3 9^2 


( 2 ) 


(3) 


and becomes identically zero for isotropic flow. The effects of rotation comes 
through the anisotropy of turbulence. In the process of turbulence modeling, if 
local isotropy of turbulence structure is assumed, the effects of rotation can be 
neglected. But, in the rotor wake, the turbulence is not isotropic. As proposed 
by Rotta, the deviation from isotropy can be written as [6] 


u . u . -ykS . . 

X j 3 xj 


(4) 


and the effects of rotation can be represented through the followxng term 


‘^1*' 2k^ 

where Cj^ is universal constant. If we introduce Rossby number (U/Lfi) in Eq. 5 
it can be also represented as follows 


(5) 


i'^k 1b —1 

*^l^~2k /k) (Rossby number) (6) 

The modeling of other terms in dissipation equation in rotating coordinate system 
can be successfully carried out much the same way as in stationary coordinate 
system. If we adopt the modeling of Eq. 1 by the Imperial College group [7,8], 
the dissipation equation in rotating coordinate system can be shown to be as 
follows 

Dc_ a.,k 3c, 

Dt ^ei k 9xj^ ^£2 1^>- ‘^e9x^1^k“£ 9x^^ 

u . I 2 

‘^1^'^k T (Rossby number) (7) 
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The effects of streamline curvature in the flow are included implicitly in the' 
momentum conservation equation. If the turbulence modeling is performed with 
sufficient lower order terms, the effects of streamline curvature can be included 
in the modeled equation. But for the simplicity of calculatxon, the turbulence 
modeling usually includes only higher order terms . Therefore further considerations 
of the effects of streamline curvature is necessary. In the modeling process, 
the turbulence energy equation is treated exactly with the exception of the diffusion 
term. The turbulence dissipation equation is thus the logical place to include 
the model for the effects of streamline curvature. 

Bradshaw [4] has deduced that it is necessary to multiply the secondary 
strain associated with streamline curvature by a factor of ten for the accounting 
of the observed behavior with an "effective viscosity" transport model. In the 
process of modeling of turbulence dissipation equation, the lower order generation 
term has been retained for the same reason. In our modeling the effects of 
streamline curvature are included in the generation term of turbulence of dissipatxon 
equation through Richardson number defined as follows. 


2Vq cosa .2. 

\ C 2 ay ^ 


( 8 ) 


Hence the final turbulence dissxpation equation which include both the rotation 
and streamline curvature effects is given by 

De _ n 4 . r R ^^i j. 1^ s 

Dt k ^^£2. k 3x^e 'Vl ^ ^l*- 2k ~ S'^ij^^k^ 


x(Rossby number) 

where Cg^, c^,, Ce 2 s ^£! are universal constants. 


(9) 


2.1.2 Turbulence closure 

The Reynolds stress term which appear in the time mean turbulent momentum 
equation should be rationally related to mean strain or another differential 
equation for Reynolds stresses should be introduced to get closure of the governing 
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equations . For the present study, k-s model which was first proposed by Harlow 
and Nakayama [9] is applied and further developed for the effects of rotation and 
streamlxne curvature. In k-s model, the Reynolds stresses are calculated via the 
effective eddy viscosity concept and can be written as follows : 


3U. 9U. 

, 1 , j, 2k, 

'^i^j 3x_,^ ~3^ij 


(- 10 ) 


J 1 

The second term in the right hand side makes the average nomal components of 
Reynolds stresses zero and turbulence pressure is included in the static pressure 
Effective eddy viscosity is determined by the local value of k & e. 


y - c — 

t p e 

where c^ is an empirical constant. The above equation implies that the k-e 


( 11 ) 


model characterizes the local state of turbulence with two parameters k and c. 
e is defined as follows with the assumption of homogeneity of turbulence structure. 


9u. 9u. 
- 1 1 


( 12 ) 


For the closure, two differential equations for the distribution of k and e are 
necessary. The semi-empirical transport equation of k is as follows 


“ = — -^1 + G - £ 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(13) 


Where G is the production of turbulence energy and G^ is empirical constant. 

In Section 2.1.1, the transport equation of e has been developed for the 
effects of rotation and streamline curvature. 


2.1.3 The numerical analysis of turbulent rotor wakes in compressors 

The rotor wakes in a ‘compressor are turbulent and develops in the presence 
of large curvature and constraints due to hub and tip walls as well as adjoining 
wakes. Rotating coordinate system should be used for the analysis of rotor 
wakes because the flow can be considered steady in this system. The above 
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mentioned three-dimensxonality as well as complexity associated with the geometry 
and the coordinate system make the analysis of rotor wakes very diffxcult. A 
scheme for numerical analysis and subsequent turbulence closure is described 
below. 

The equations governing three-dimensional turbulent flow are elliptic and 
large computer storage and calculation is generally needed for the solution. The 
time dependent method is theoretically complete, but if we consider the present 
computer technology, the method is not a practical choice. The concept of 
parabolic-elliptic Navier-Stokes equation has been successfully applied to some 
simple three dimensional flows. The partially parabolic and partially elliptic 
nature of equations renders the possibility of marching in downstream, which 
means the problem can be handled much like two-dimensional flow. The above 
concept can be applied to the flow problems which have one dominant flow direction. 
Jets and wakes flows have one domxnant flow directxon and can be successfully 
solved with the concept of parabolic elliptic Navxer-Stokes equation. The 
numerical scheme of rotor wakes with the concept of parabolic elliptic Navier- 
Stokes equation are described in Appendix II. Alternating direction implicit 
method is applied to remove the restrictions on the mesh size in marching direction 
without losing convergence of solution. The alternating-direction implxcit method 
applied in this study is the outgrowth of the method of Peaceman and Rachford 
[ 10 ]. 

Like cascade problems, the periodic boundary condition arises in the rotor 
wake problem. It is necessary to solve only one wake region extending from the 
mid passage of one blade to the mid passage of other blade because of the 
periodicity in the flow. This periodic boundary condition is successfully 
treated in alternating direction implicit method in y- implicit step (Appendix II]. 
It is desirable to adopt the simplest coordinate system in which the momentum 
and energy equations are expressed. As mentioned earlier the special coordinate 
system should be rotating with rotor blade. 
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The choice of special frame of the reference depends on the shape of the 

boundary surface and the convenience of calculation. But for the present study, 

the marching direction of parabolic elliptic equation should be one coordinate 

direction. The boundary surfaces of present problem are composed of two concentric 

cylindrical surfaces (if we assume constant annulus area for the compressor) and 

two planar periodic bomdary surfaces. One apparent choice is cylindrical polar 

coordinate system, but it can not be adopted because marching direction does not 

coincide with one of the coordinate axis. Ro simple orthogonal coordinate system 

can provide a simple representation of the boundary surfaces (which means boundary 

surfaces coincide with the coordinate surface, for example, x^ = constant) , In 

view of this, the best choice is Cartesian because it does not include any terms 

arising from the curvature of the surfaces of reference. The overall description 

is included in Appendix II in detail. OBIGINAL PAGE IS 

OP POOR QUALITY 

2 . 2 Modified Momentum Integral Technique 

This analysis* is a continuation and slight modification to that given in 
Reference 11. 

In Reference 12, the static pressure was ass\amed to be constant across the 
rotor wake. From measurements taken in the wake of heavily and lightly loaded 
rotors a static pressure variation was found to exist across the rotor wake. 

These results indicate that the assumption of constant static pressure across 
the rotor wake is not valid. The following analysis modifies that given in 
Reference 11 to include static pressure variation. 

Also included in this analysis are the effects of turbulence, blade loading, 
and viscosity on the rotor wake. Each of these effects are important in controlling 
the characteristics of the rotor wake. 

*Even though Equations 14 through 43 given in this report are a slightly modified 
version of the analysis presented in Reference 11, the entire analysis is 
repeated for completeness. The analysis subsequent to Equation 43 is new. 
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The equations of motion in cylindrical coordinates are given below for r, 0, 
and z, respectively (see Figure 1 for the coordinate system) 




_ 9(uV) 

p 3r 3 r r- r ” r30 


8(u'w') 

9z 


V 9V , „ 3V 9V UV _ 1 9P 1 9(v'^) 2(u'v’) 9(u'v') 

r96+^9r + ^9z^ r ^ ~ ~ ~ x ~ ? ^ 


(14) 


9(vV») 

9z 


IM+U — +W-^=-i-^ - _ 3(u*w') _ (uV’) _ 1 9(v'w') 

r 90 ^ 3r ^ 9z " p 9z ” 3z 9r ~ r r 90 

and continuity. 


9V ■ ^ , 
r90 9z 


9U 

9r 



(15) 

(16) 

(17) 


2.2.1 Assumptions and simplification of the equations of motion 

The flow is assumed to be steady in the relative frame of reference and to 
be turbulent with negligable laminar stresses. 

To integrate Equations 14-17 across the wake, the following assumptions are 

made, 

1. The velocities and are large compared to the velocity defects u 
and w. The wake edge radial velocity, U^, is zero (no radial flows exist outside 
the wake). In equation form this assumption is, 

U = u(0,z) (U^ = 0) (18) 

V = V^(r) - v(0,z) (V^ » v) (19) 

W = W^(r) - w(0,z) (W^ » w) (20) 

2. The static pressure varies as, 

p = p(r,0,z) 



3. Similarity exists in a rotor wake. 

Substituting Equations 18-20 into Equations 14-17 gives. 


9 


o 

r 


IH 

30 


+ W 


3u 

32 


- -(V 
r o 


+ Qv) - 


-(-2vV - 2vfJr) 
r o 


2 2 2 

1 lE _ 9^' _ + ill 

p 3r 9r r r 


3(u'v') _ 3CuVM 
r30 3z 


( 21 ) 


r 30 


u 


3V 
o 

3r 


W 

o 


3v 

3z 


+ 2fiu + 



1 3p _ 1 3(v*^) _ 2(u'v') _ 3(u'v') 
p r39 ~ r 30 ~ r ~ 3r 


3(vV») 

3z 


Z° ^ ^ ^ _ 3 (w' _ 3 (u'w* ) _ (u*w* ) _ _1 3 (v*w* ) 

r 30 3r o 3z p 3z 3z 3r r r 36 


( 22 ) 

(23) 


2. Jx _ lx J. Ji = 

r 30 3z r 


(24) 


2.2.2 Momentum integral equations 

Integration of inertia terms are given in Reference 12. The integration of 
the pressure gradient, normal intensity and shear stress terms are given in this 
section. 

The static pressure terms in Equations 21-23 are integrated across the wake 
in the following manner. For Equation 22, 


^®o 1 3p j _ 1- 

r 30 r^^e 

c 




The integration of the pressure term in Equation 21 is, 


r o Ip 

0^ 3r 
c 


d0 



36 30 

c o 

p — _ — p 

c 3r e 3r 


It is known that, 


r(0 - 0 ) = 5 

c o 

Differentiating the above expression. 


■jjCre^ - r8^> = 


12 

3r 


(25) 


(26) 


(27) 


ORIGINAL 

OF POOR (Quality 


(28) 
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r99 „ r30 

O _ 1^.4. Sl 

3r r ~ 3r 3r 


(29) 


Using Equation 29, Equation 26 becomes, 

0 


6 


f 




36 


9, de (p^ - p^)-^ - p^(A _ 

r 


(30) 


Integrating the static pressure term in Equation 23 across the wake. 
6 . - p 30 36 

‘iTN ^ rs /- 


o l£ HO = iE. 4- 

Q' 3z ^ 3z 3z 

c 


J 


~ p 


e 3z 


(31) 


where 


36 


3z 


and 


c ^ jL .V ^ 1 ,^0 ~ ^c 

rSj'^c r^W - w 
o c 


If = i^-^c - -^o) 


30 30 . 

o c 1 35 


•) 


(32) 


3z 9z r 3z 

Substituting Equations 32 and 34 into Equation 31 


(33) 

(34) 


6^ 5 p V - V 

/ o i£ dz = -^ + ? 

6 3z 3z ^ r ^ W„- w 

C o c 


1 V - V 

.) _ p [A(_2 ] 

^ e^r^W„- w ^ r Sz-* 


1 35. 


(35) 


Integration across the wake of the turbulence intensity terms in Equations 21- 
23 is given as, 


^ .1 , 3u*~^ 

r o 3r 


— + — )dn 
r r 

(36) 

30 

(37) 

,2 


3z 

(38) 


T / (- 


Integration of any quantiy q across the wake is given as (Reference 12) , 


d0 = -a 

6-^ 36 

c 


(39) 
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6 


d - o 


30 


0 


; O de = ^ - q(z,0)d9 + — - q. 


30 


3z 


c c 

If similarity in the rotor wake exists 

q = q (z)q(n) 


‘c 3z 3z 


(40) 


(41) 


Using Equation 41, Equation 40 may be written as, 

0 
I 

0 


/ ° 3(q(z,0)) ^ ^1 

c 3z rdz o 


q V - V 

. r q(q)dq + — (-^ £•) 

rdz o r W - w 


(42) 


From Equation 40, it is also known that, 

0 6q ^ 

° qd0 = / q(n)dn (43) 

The shear stress terms in Equation 21 are integrated across the wake as 
follows. First, it is known that 



3(u'w') 

r30 


d0 



3(u^v’) 

r30 


d0 = 0 


(44) 


c c 

since correlations are zero outside the wake as well as at the centerline of the 
wake. The remaining shear stress term in Equation 21 is evaluated using the eddy 
vxscosity concept. 


U » - 


Neglecting radial gradients gives, 


I I /^n\ 


Hence, it is known that. 




. 0 
d . 0 


0 


dz - -Pt ^ ° If 

C C 


Using Equation 40 and Equation 47, 


r^o 8 . , ,, ,, 3rd JO ^ 

c c 


(45) 


(46) 


(47) 


( 48 ) 
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^ TT J U V - V 

9 fH d. / f. -y , HX/ o Cy -i 

0 z r dz^ r ~ vr ^ 


(49) 


Hd/p V Akip 

■ -'^T ¥ ttC'V - — — '■ 

dz 


y„ j V u V u V w^u 
TdrOm cm. ocm 


r dz"- W 


W 


+ 


] (50) 




Neglecting second order terms. Equation 50 becomes, 


9 « lirr, j2 V d(u_) 

-/ ° + iT TT-1 

c dz ‘ o 


(51) 


H is defined in Equation 58. Using a similar procedure the shear stress term 


in the tangential momentum equation can be proved to be, 

^WTae - - ^[G ^(av^) + ^ ^1 


(52) 


G is defined in Equation 58. The radial gradients of shear stress are asstimed to 
be negligible. 

Using Equations 39, 42, and 43, the remaining terms in Equations 21—24 are 
integrated across the wake (Reference 12) . The radial, tangential, and axial 
momentum equations are now written as. 


g <i(% 4 ) (V^ + nr) 5 S 26 G(V^ + nr)v^S 

r ^ “di O T “2 


r W 


W r 
o 


1 JE A _ _L s !!s.(p - p ) + !!a(X - A1 
P ar 2 „2 “ 3r ' c „2 ,2 rdr 

Wo PW„ ^ 


,) +i.A dt-lSL- 

’ * T 3r 

o 


,2 


f2 ,2 

n . V Tj 

— — + — —-Jdn 


Sy™ .2 V du 

r~[E -~( 6 u ) + J] 
W r ,2 m W dz 

o dz . o 


(53) 


5u H dV V . 

r. m r o . o . GS d , 

S ■ ; — [-n — + — + 2£2] 3 — (v 5) = 

W r '• dr r ■' r dz c 


— — S (P - P ’ ) + 

„2 ^''^e ^c^ ^ 7 o 

WoP W^r 


2 

,1^ av’ ^ , 

f (- 


r96 


Sy„ ,2 V dv 

^tG 


( 54 ) 
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dr W^r - r dz ^ " ^,2^ 3z pr ' c dz 


2 

, <5S ^Ir Sw’ TJ 

+ -:t- f [ — rr"idn 


W^r o 
o 


3z 


and continuity becomes 

Sv 


d(6w ) S5u H Sw V 
_c _ ^ c^ m c o ^ ^ 

r r dz 2 rW 


where w , v » and u are nondimensional and from Reference 12, 
c c m 


G = /^g(n)dn , 


H = /\(n)dn , — y Y = h(n) 




(55) 


(56) 


(57) 

(58) 

(59) 


2,2,3 A relationship for pressure in Equations 53-55 

The static pressures and p are now replaced by knoxm quantities and the four 

unknoTTOS v , u , w , and 6. This is done by using the n momentum equation (s, n 
c m e 

being the streamwise and principal normal directions respectively - Figure 1) of 
Reference 13, 


U 


3U 
n 

't 3r 


+ U 


3U 


n 


n 


8n 


U 


3U 
n 

3s 


+ 2S^U cos 3 - 


U" U U . 

s . r n 2, 

1 cos fc 


R 


1 Bf* r 3 ^ r t N , , 2^^^^r n , 

X xr— (u'u ) + (1 + cos 3) + 

p 3n 3r^ r n^ ^ r 


3 ,2^3- 

— u' + — 

an n 


. (u'u’) 
3s n s 


(60) 


where 


r, 1 2 

P — y r 


(61) 


Non-dimens ionalizing Equation 60 with S and U„ and neglecting shear stress terms 

^o 

results in the following equation. 
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Ur \ Ug *S(U^/Us^) 26 (fir) Ur cos 5 UrU^ 2. 

UZ 6(r/6) UT -S(n/5) iT" 6(s/5) ^ 


Uso’^ Ug^(r/6) 


U^/Us 

s 

R /6 

c 


1 9P^ 


p 9(n/6) ^2 


3 (n/6) 


(62) 


An order of magnitude analysis is performed on Equation 62 using the following 
assumptions where e « 1, 


n V 6 

Un/Uso ^ 

s V 6 

U^/Usq V 1 

6/r V e 

u'2^/u2 V 
n so 

fir/Us 1 

R V S 

‘’o 

c 

U /Us '^1 
s ®o 

S/S 'V' e 


in the near wake region where the pressure gradient in the rotor wake is formed. 


The radius of curvature, (Figure 1) is choosen to be of the order of blade 
spacing on the basis of the measured wake center locations shoxm in Figures 2 
and 3. 

The order of magnitude analysis gives. 




\ 3 cn„/u.„) 2 

V e 


Ug^ 3 (n/6) 


aCuZUs,) 


Ug^ 3(s/6) 


o> g 


26 (fir) U cos 6 

— 

^o 


V e 


(63)> 


(64) 

(65) 


( 66 ) 


2,2 

cos 6 ^> £ 


Uso(r/6) 


(67) 
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U^/Ug 

S °Q 

R /6 

c 


l\j £ 


( 68 ) 


3(u;^/Ug\) 


e 


3(n/5) 

Retaining terms that axe of order e In Equation 62 gives 


2(i2r)u^ cos g U^/u|^ ^ ^ 




3s 


Ug r 
So 


R 


P 3n J.2 

Us, 


3n 


(69) 


(70) 


The first term in Equation 70 can be written as, 

0, acn^Hs,) 3((B„/up (n^/Hs^)) 


u. 


3s 


Us 


3s 


(71) 


Since, 


TJ /TJ= = tan a 
n ®o 


where a is the deviation from the TE outlet flow angle. Equation 71 becomes, 


U 3(U„/Ug J 


n 


Uc 


3s 


U , U 

S V r i_ St 

uT 


“s “s 2 3.,, 3(D^/t),„) 

■jt-toT “-sS- “ 5S — J 

==o ‘’O 


U , U 3(U /Uc. ) 

S 1 , S S Sq 

„2 R “ Ug 3s 

Uso ^ 


(72) 


where 

3a _ ^ 
3s R 

c 


and 


2 , 
sec a = 1 

tan a a 

with a small. Substituting Equation 72 and Equation 61 into Equation 70, 
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U 9(U /Us„) 2(S^r)U cos g 


3 S 


IDP 1 

P 9n ^2 an 


The integration across the rotor wake in the n direction (normal to streaim^ise) 
would result in a change in axial location between the wake center and the wake 
edge. In order that the integration across the wake is performed at one axial 
location, Equation 73 is transformed from the s, n, r system to r, 9, z. 

Expressing the strearawise velocity as 

2 2 2 
U = V + W 
s 


^ - cos g , — = sxn g 


ORIGINAL PAGE IS 
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Equation 73 becomes, 




2 2 

V + W 
o o 


o o 


o o 


2 2 

V + W 
o o 


1 cos g e c „ 9 / n . . ^ 9 ;r-N 

" p r ^ ^ ^ 92 ^ 

O ■ o ^®o ° 

Expressing radial, tangential, and axial velocities in terms of velocity defects 
u, V, and w respectively. Equation (74) can be written as, 

1 ar cos 6r , rj 3v, ^®o,3w^ „ 3w^.„^ la . . 

2 ^2 ^ le^ '-37 - 2 ;r-V <1F + 

0.0 o "o 

6 2 0 
9v. j . o.9w , 9Wv , - 2(Qr)r rO,. 

+ ^o + 0^ + 2 6^ = 

c V + W c 

o o 


2 2 

1 - sen 8 '75) 


V + w 
o o 


c u: 


Using the integration technique given in Equations 39-43, Equation 75 becomes, 
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d(v^6) d(v 6) V - V d(w^5) d(w S) 

o [ — ;3 G- + V — ~ — G + (v^ + V V ) (-^ -) + — F + W ^ — F 

dz 1 o dz ^ '■''c o - w/ dz 1 o dz 


o dz 


o c 


P V - V W V 

+ (w + w W ) ] + ~ tan 6[v + v + + w ] 

C COW_^ — w 2 cWc C c 

0 c o 


,.2 


(W^ + V^) 9 ^ u'^ 0 , u- 

-t 2w^(.r)du^H -t ° A(J_)ae H- ° A(-f-)de 


u= 


Uo 


= - p<"e; - 


(76) 


where. 


V 


S S^<Tl)dr! , g(n) , , 

J- o 


5’-, - f^(n)dri , f(n) == 

j- r\ 


W 


w^(z) 


(77) 

(78) 


Non-dimensionalizing Equation 76 using velocity W gives. 


d(v^5) 


V. d(v 5) 




_ +irv:iv .:ir(w„6)F 


V - V d(x/6) 


dz 1 W dz 
o 


c W V - w 

o o c 


dz 1 


F, + dz^”c^ 


■c ■ "c^ 3 ^ 1 e [v^ + + w? + wj 


+ (w_ + w^) (- 


o c 


c W c c c' 
o 


2 2 2 2 

+ 2-%^ u «H + (1 + -4)./ ° — (-4-)ae + (1 + — ) /° — (-^-)d6 

o “o o ®o 


-- f-(p - p ) 
,,2 e c 
W p 
o 


(79) 


Note, v^, w^, and u^ are now normalized by W_^ in Equation 79. 


2.2.4 Solution of the momentum integral equation 

Equation 79 represents the variation of static pressure in the rotor 
wake as being controlled by Corriolis force, centrifugal force, and normal gradients 
of turbulence intensity. Prom the measured curvature of the rotor wake velocity 
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centerline shown in Figures 2 and 3, the radius of curvature is assumed negligible. 
This is in accordance with far wake assumptions previously made in this analysis. 


Thus, the streamline angle (a) at a given radius is small and Equation 79 becomes. 


25(S!r)uH e , u'^ 7^ e - 

O W C TJa W c U«. 

O o O ®0 


(p - p ) 
e c 


Substation of Equation 80 into Equations 53, 54, and 55 gives a set of four 
equations (continuity and r, 9, and z momentum) which are employed in solving for 
the unknoxTOS 5, u , v , and w . All intensity terms are substituted into the final 
equations using known quantities from the experimental phase of this investigation. 

A fourth order Runge-Kutta program is to be used to solve the four ordinary 
differential equations given above. Present research efforts include working tox'/ards 
completion of this programming. 
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3. EXPERIMENTAL DATA 

3.1 Experxmental Data From the AFRF Facility 

Rotor wake data has been taken with the twelve bladed un-canibered rotor 
described in Reference 13 and with the recently built nine bladed cambered rotor 
installed- Stationary data was taken with a three sensor hot wire for both the 
nine bladed and tweleve bladed rotors so that mean velocity, turbulence intensity, 
and turbulent stress characteristics of the rotor wake could be determined. 

Rotating data using a three sensor hot wire probe was taken with the twelve bladed 
rotor installed so that the desired rotor wake characteristics listed above would 
be known. 

3.1.1 Rotating probe experimental program 

As described in Reference 1, the objective of this experimental program is 
to measure the rotor wakes in the relative or rotating frame of reference. Table 1 
lists the measurement stations and the operating conditions used to determine the 
effect of blade loading (varying incidence and radial position) and to give the 
decay characteristics of the rotor wake. In Reference 1 the results from a 
preliminary data reduction were given. A finalized data reduction has been completed 
which included the variation of due to temperature changes, wire aging, and 
other changes in the measuring system which affected the calibration characteristics 
of the hot wire. The results from the final data reduction at r/r^ = 0.720, i = 

10°, and x/c = 0.021 and 0.042 are presented here. Data at other locations will 
be presented in the final report. 

Static pressure variations in the rotor wake were measured for the twelve 
bladed rotor. The probe used was a special type which is insensitive to changes 
in the flow direction (Reference 1) . A preliminary data reduction of the static 
pressure data has been completed and some data is presented here. 
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Table 1 

Experimental Variables for Rotating Hot Wire and 
Static Pressure Measurements 

Incidence 0®, 10° 

Radius (r/r^) 0.720 and 0.860 

Axial distance from trailing edge 1/8", 1/4", 1/2", 1", and 1 5/8" 

Axial distance x/c (non-dimensionalized 0.021, 0.042, 0.083, 0.167, and 

w.r.t. the local chord) 0.271 

The output from the experimental program included the variation of three 
mean velocities, three r.m. s. values of fluctuating velocities, three correlations, 
and static pressures across the rotor wake at each wake location. In addition, 
a spectrum analyzer was utilized to derive the spectrum of fluctiiiating velocities 
from each wire. 

This set of data provides properties of about sixteen wakes, each were defined 
by 20 to 30 points. 

3. 1.1.1 Interpretation of the results from three-sensor hot-wire final data reduction 
Typical results for r/r^ = 0.72 are discussed below. This data is graphed 
xn Figure 4 through 24. The abscissa shows the tangential distance in degrees. 

Spacxng of the blade passage width is 30°. 

(a) Streamwise component of the mean velocity, US 

Variation of US/USO downstream of the rotor wake. Figures 4 and 5, indicates 
the decay characteristics of the rotor wake. The velocity defect at the wake 
centerline (USC)^/US0 is about 0.45 at x/c = 0.021 and decays to 0.56 at x/c = 0.042. 

The total velocity defect (Qd)^/Q^ at the centerline of the wake is plotted 
with respect to distance from the blade trailing edge in the streamwise direction 
in Figure 6 and is compared with other data from this experimental program. Also 
presented is other Penn State data, data due to Ufer [14], Hanson [15], isolated 
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airfoil data from Preston and Sweeting [16] , and flat plate wake data from 
Chevray and Kovasznay [17]. In the near wake region decay of the rotor wake is 
seen to be faster than the isolated airfoil. In the far wake region the rotor, 
isolated airfoil, and flat plate wakes have all decayed to the same level. 

Figures 4 and 5 also illustrate the asymmetry of the rotor wake. This 
indicates that the blade suction surface boundary layer is thicker than the 
pressure surface bomdary layer. The asymmetry of the rotor wake is not evident 
at large x/c locations indicating rapid spreading of the wake and mixing with 
the free stream. 

(b) Normal component of mean velocity, UN 

Variation of normal velocity at both x/c = 0.021 and x/c == 0.042 indicate 
a deflection of the flow towards the pressure surface in the wake. Figures 7 and 
8. The deflection is seen to decay from x/c = 0.021 to x/c = 0.042. This decay 
results from flow in the wake turning to match that in the free stream. 

(c) Radial component of mean velocity, UR 

The behavior of the radial velocity in the wake is shown in Figures 9 and 10. 
The radial velocity at the wake centerline is relatively large at x/c = 0.021, 
and decays by 50 percent by x/c = 0.042. Decay of the radial velocity is much 
slower downstream of x/c = 0.042. 

(d) Turbulent intensity, FWS, FWN, FWR 

Variation of FWS, FWN, and FWR across the rotor wake at x/c = 0.021 and 
x/c = 0.042 as shown in Figures 11 and 12, Figures 13 and 14, and Figures 15 and 
16, respectively. Intensity in the r direction is seen to be larger than the 
intensities in both the s and n directions . Figures 15 and 16 show the decay of 
FI-TR at the wake centerline by 50% over that measured in the free stream. Decay 
at the wake centerline is slight between x/c = 0.021 and x/c = 0.047 for both 
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FWS and FHN. At x/c = 0.167, FWS, FWN, and FWR have all decayed to close to that 
in the free stream. 

This final data reduction indicates large gradients of, 

9 (FWR) 9 (FWR) 

9z ’ 90 

with much slower decays of FWS and FWN in the near wake. To correctly predict 
the flow in this region, the terms' containing large gradients of FITR must be 
retained. 

(e) Turbulent stress, B(4), B(5), and B(6) 

Stresses in the s, n, r directions are shown in Figures 17 and 18, Figures 
19 and 20, and Figures 21 and 22, respectively. Turbulent stress B(4) is seen 
to be negative on pressure surface side of the wake centerline, positive on the 
suction surface side, and zero at the wake centerline. Decay of B(4) from x/c = 
0.021 to x/c = 0.042 is large. Beyond x/c = 0.042 decay of turbulent stress B(4) 
is much slower. Turbulent stress B(5) is small compared to the levels measured 
for B(4) and B(6) . Turbulent stress B(6) is 0.015 at x/c - 0.021 and decays at 
.05 of the peak seen on the pressure surface side in the rotor wake at x/c = 0.042. 
Decay is much slower do^stream of x/c = 0.042. 

(f) Resultant stress, EESSTR 

Variation of resultant stress across the rotor wake is shorn in Figures 23 
and 24 for x/c = 0.021 and x/c => 0.042. Resultant stress is highest on the 
pressure surface side of the wake centerline. This peak decays from a resultant 
stress level of 0.017 at x/c = 0.021 to a level of 0.005 at x/c = 0.042. At x/c = 
0.271, RESSTR has decayed to close' to its free stream value. 

3.1.1. 2 Interpretation of the results from static pressure preliminary data 
reduction 

Variation of static pressure across the rotor wake is shoxra in Figure 25 
through 27 at r/r^. = 0.721 for x/c = 0.021, 0.083, and 0.271. Local static 
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pressure in the wake is non- dimens zonal ized with wake edge static pressure, P , on 
the pressure surface side of the wake. Large gradients of P/P^ is 'clearly seen 
at all axial locations. This indicates that when theoretically predicting the 
rotor wake that terms such as, 

9P 

90 ’ 9z 

must be retained. However, the variation in static pressure shown in Figure 25 
through 27 is not clearly understood. At a larger radial position (r/rj- = 0,860), 
the static pressure variation shoxTO in Figure 28 was measured at x/c = 0.021. 
Comparison of Figures 25 and 28 does not show a similar variation of P/P^ across 
the wake as would be expected. For this reason, the results shown here are only 
preliminary as the data reduction procedure used in calculating P is presently 
being checked. Therefore, the only reliable conclusion that can be drawn from 
this preliminary data reduction is that gradients of static pressure do exit in 
the rotor wake. 

3.1.2 Stationary probe experimental program 

Using a three-sensor hot wire probe, stationary data was recorded using the 
technique described in Reference 13. Measurement locations and the operating 
conditions used in the measurement of the wake of a-12 bladed uncambered rotor is 
described in Referenced. With a 9 bladed cambered rotor installed, stationary data was 
taken at three operation conditions in both the near and far wake regions. Axial 
distance was varied from x/c = 0.042 to about 1.2 chords downstream at ten radial 
positions from hub to tip = .465, .488, .512, .535, .628, .721, .814, .916, 

.940, and .963). This data will indicate the effect of end wall flows and the 
hub wall boundary layer on the rotor wake. To determine the centrifugal force 
effect on the rotor wake, *data was recorded at' design conditions with two rotor 
RPM's (RPM = 1753 and 1010). To indicate the effect of loading, the data was also 
taken at an additional off-design condition at each measurement location. 
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The nine bladed data described here and the previously taken twelve b laded 
data will be reduced when the data processing computer programming has been 
completed 

3.2 Experxmental Program at the Axial Flow Compressor Facility 

Rotor wake measurements to be carried out at the axial flow compressor facility 
consist of both the rotating and the stationary hot wire measurements along with 
a static/stagnation pressure survey across the wake. The purpose of this 
investigation is to study the following: 

1. Wake development and decay as it travels downstream. 

2. Nature of xcake along the radial direction. 

3. .Nature of x^ake in the hub and annulus xjall boundary layers. 

4. Effect of blade loading on the wake. 

5. Effect of IGV wakes on the rotor xo-ake. 

6. Effect of inlet distortions on the x^rake. 

7. Comparison of measurements in stationary and rotating frame of reference. 

For the pressure survey a special type of static/stagnation probe is being 

fabricated. It is designed to be insensitive to direction changes. 

Presently, the rotating hot x^ire measurements are being carried out," Using 
the traverse gear, which x^as modified to suit the present requirements, the three 
sensor hot xfire is being traversed across the wake. A sketch of the modified 
tranverse gear is shoxm in Figure 29. The output which is in the form of three 
D.C. voltages, three r.m. s. voltages and three correlating voltages are being 
processed. Some preliminary results of the rotor xjake, measured at radius = 

9.448" (R/R^ = .45) and 1.75 inches axial distance doxmstream of the trailing 
edge (corresponds to 0.3 times chord length) at the design operating condition 
= 0.54) are shoxm in Figures 30 and 31. The xsridth of the passage is 18°. The 
streanEvrise velocity US, shoxTO. in Figure 30, shows a rather slower decay of the 
wake in this configuration of heavily loaded rotor. The radial velocities, shoxm 
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in Figure 31, are also higher than those reported in previous cases. The 
turbulence intensities (not shoxra here) are also found to be high. 
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4, SYNTHESIS AND COERELATION OF ROTOR WAKE DATA 

4.1 Similarity rule in the rotor wake 

There is sufficient experimental evidence about similarity of mean velocity 
profiles in successive sections of wake flows in the doT-rastream direction. The 
free— turbulent— flow regions are relatively narrow, with a main— flow velocity 
much greater than the transverse component. The spacial changes in the main- 
flow direction are much smaller than the corresponding changes in the transverse 
direction. Furthermore, the turbulence flow pattern moving do’cmstream is strongly 
dependent on its history. These facts make similarity probable in wake flows. 

The similarity rule is examined for rotor wake flows. The data due to Raj and 
Lakshminarayana [13] and Schmidt and Okiishi [18] as well as those present in 
Section 3 are examined, 

4.1.1 Data of Raj and Lakshminarayana [13] 

The maximum mean velocity difference is used as a velocity scale to make all 
the velocities non-dimensional. The so-called "half— value" distance which is the 
distance from the axis of symmetry to the locatxon where the mean velocity defect 
is half the maximum value, xs used as a length scale. 

Similarity in axial velocity profile was discussed in Reference 12 and will 
not be repeated here. 

Due to imbalance between centrifugal force and static pressure in the radial 
direction, the radial velocity component is formed. In the blade boundary layer, 
the radial velocity profile is determined by the above mentioned imbalance of 
forces acting on the fluid particle. The experimental data of rotor wake shows 
that there are faro dominant radial velocity scales. One is dominant in the center 
region and the other is dominant in the outer region of the wake. Figure 32 
shows general profile of radial velocity and notations used. In the center 
region (0 _< |y| < jL^j or Ih^j) the velocity is non-dimensionalized with (Auj-^)p or 
(Au;i-^)s. In the outer region the velocity is non-dimensionalized by (Auj.)p or 
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(Au ) • The length L or L is used for non-dimens ionali zing length scale. 

P s P s 

Figure 33 and Figure 34 show typical radial velocity distribution with the above 
mentioned velocity and length scales. In the center region, the profile is 
nearly linear. In the outer region, the velocity is very small, and hence the 
large scatter in the data. 

4.1.2 Data of Schmidt and Okiishi [18] 

The velocity profiles with the same corresponding velocity and length scales 
are shown in Figure 35 and Figure 36 for the Iowa State data. Better agreement 
with theoretical distribution in axial velocity component is shorn. Good agreement 
with theoretical Gaussian distribution is also shown for the absolute tangential 
velocity profile plotted in Figure 37. 

With regard to the radial velocity component, plotted in Figure 38, the 
new velocity and length scale defined earlier are used successfully. The bell- 
shaped velocity profile curves are found in axial and tangential velocity components. 
But for radial velocity components, different length and velocity scales are 
necessary and the velocity profile is nearly linear near the center region- 

4.1.3 Rotating probe data from the AFRF facility 

The data reported in Section 3 of this report is plotted in Figures 39 through 

50 in a non-dimensional form using the characteristic length scale (L and L ) 

p s 

and velocity (w^, u^ and u^) discussed earlier. These results are from a twelve 
bladed rotor operating at 0° and 10“ incidence and taken with a rotating triaxial 
probe. 

The axial velocity profile, shown plotted in Figure 39 and 40, seem to follow 
the Gauss’ function (e tangential velocity defect, shown in 

Figure 41 and 42, show the same trend. Even the radial velocity seem to folloxvf 
this trend (Figure 43) . The nature of thse radial velocxty profiles are different 
from those reported earlier (Figures 33, 34, and 38). Further work x^ith regard 
to the existance of similarity in radial velocity is presently underxj'ay. 
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The radial, axxal and tangential component of turbulence intensity profiles 
at various incidence, radial and axial locations are shoxim in Figures 44 through 
49 . The intensities are plotted as 


,2 

,2 u’ 

U B 

VS n or n 

-;2 ^ p " 


The intensities are normalized with respect to the extrapolated peak intensity 


near the center and the corresponding value in the free stream is substacted from 
the local value. All the data seem to follow the Gauss' function 


4.2 Fourier Decomposition of the Rotor Wake 

The rotor wake can be conveniently represented in terms of a Fourier Series. 
The velocity at a position (r, 0, z) can be represented as; 


11 ^ -v* 


(81) 


where 


u. 


Ujj = normalized velocity defect = 


D 


n 


V 


u - u 
max 

u “ u . 
max min 


8^ = distance between wake centerline and the location where the 
^ 1 

= -2 (x^ake width) ’ 

u = maximum velocity in the wake 
max •' 

u , = centerline velocity in the wake 

min ■’ 

u = velocity at any point 

In the present study, the wake data from the axial flow research fan reported 
in Reference 13 is analyzed by Fourier decomposition. The method is based on the 
recursive technique of Ralston [19]. In the analysis, the velocity defect was 
normalized by the wake centerline defect and all the angles by the angle at the 
half wave x^idth in Equation 81. represents half the average of the wake velocity 
in the given interval. 
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Now 


n 


fn 

2ir o Uj) 


cos 


n9 

e. 


de 


where e. = Nemnann's factor 1 for n = 0; = 2 for n > 0) 
n 

2tt = interval of integration (in the analysis only the wake portion is 
considered) . 

1 1 2tt 

A = -^ / 

o 2 it o T 

= average over the interval 

Similarly and are calculated. As the coefficients represent the normalized 
velocity defect, it should stay constant if similarity exists. Figures 50 to 56 
shows the coefficients plotted against distance from the trailing edge normalized 
by the blade chord. It is obvious that there is considerable scatter in the 
coefficients indicating the scatter in the data. Figure 56 shows the scatter 
of the coefficients for the first four harmonics, which were found to be the 
dominant ones. The values of the amplitude of the sine component of the 
series was found to be extremely small indicating the axisymmetric nature of the 
wake, where all the distances were normalized by the respective length scale and 
the local velocity defect by the centerline defects. The Fourier series matches 
quite well with the data points, when only the first four harmonics were considered. 
Figures 57 and 59 shows how closely the two curves match. The correlation 
coefficient for the entire wake was found to be around 0.9988. 
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APPENDIX I 
Nomenclature 


B(4) ,B(5) ,B(6) =5 r.m.s. turbulent stress in the streamwise, normal, and 

radial directions respectively (normalized by local resultant 
dynamic head) 

c = Chord length 


= Lift coefficient 

c ,c ,c -,c „,c = Universal constants 

1 e el £Z c 

F,G,H = See Equations 57-59 


^ 1 ’% 


= See Equations 77 and 78 


F^S,FWN,EWR 



i 

k 


m 


=» r.m.s. turbulent intensity in the streamwise, normal, and 
radial directions, respectively (normalized w.r.t. the local 
resultant velocity) 

= Metric tensor 

= Incidence 

1 ~2 ~Z ~2 

~ Kinetic turbulent energy ^ 

= Passage-mean mass flow rate 


L ,L 
P s 


- Wake wxdth at half depth on the pressure and suction side 
respectively 


P 


= Static pressure 


P 

P' 

P" 

p* 

Q 

Qd 

RESSTR 


= Wake center static pressure 
= Wake edge static pressure 

= Nean static pressure at y,r plane (in momentum integral 
analysis thxs is wake average pressure) 

= Fluctuating component of static pressure 

= Correction part of static pressure 

= Reduced pressure 

= Total relatxve mean velocity 

= Defect in total relative mean velocity 

= Total resultant stress -/(B(4)]^ + [B(5)]^ + [3(6)]^) 

= Local radius of curvature of streamwise flow in wake. 
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r,0,z 


S 


= Rotating cylindrical coordinate system (radial, tangential, 
and axial directions, respectively, z = 0 is trailing edge. 
Figure 1) 

= Blade spacing 


s,n,r 



u,v,w. 


= Streaun^ise, normal, and radial directions, respectively 
(Figure 1) 


^ 8u. 3u, 

1 i i 

= Strain rate fluctuations + "sx ^ ^ 

= Relative radial, tangential, and axial mean velocities, 
respectively (Figure 1) 


US, UN, UR 


= Relative streamwise, normal, and radial mean velocities, 
respectively 


USO 


= Relative streamwise mean velocity in free stream 


U. 


U , U 
n’ r 


U. 

1 


u,v,w 






= Relative streamwise normal and radial mean velocities, 
respectively (Figure 1) 

= Component of mean velocity 

= Defect in radial, relative tangentail and axial velocities, 
respectively, in the wake 

= Components of turbulent velocity 

= Value of u, v, and w at node coordinate (j,k) 


u 


m 


= Maximum radial velocity away from the centerline normalized 

by W 
^ o 


u’,v',w' = Turbulence velocities in the radial, tangential, and axial 

directions, respectively (Figure 1) 

u' v' ,u'w' ,v'w' = Turbulence velocity correlation 

u*,v*,w* = Intermediate value of u, v, p during finite different 

calculation 

u',u^,u^ = Turbulent velocities in the streamwise, normal and radial 

directions, respectively (Figure 1) 

V = Relative tangential and axial velocity, respectively, outside 

of the wake 


V 


6 


^c’ 


w 


= Streamline velocity in the definition of Richardson number 

= Defect in tangential (relative) and axial velocity, 
respectively, at the wake centerline normalized by 

= Axial distance 


X 
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Ax,Ay, Az 

Q 

a 


e 

6 


5. 

£ 


e . . 

13 


e 


ijk 


e 

c 



Ti >n 

's 'p 
X 

Urp 


P 


= Mesh size in x, y and r direction, respectively 
= Angular velocity of the rotor 

= Angle between streamwise and blade stagger directxons 
(Figure 1) 

= Angle between axial and streamtirise directions (Figure 1) 

= Semi wake thickness [r(0 -0)1 

c o 

= Kronecker-delta 
= Turbulence dissipation 
= Dissipation tensor 
== Alternating tensor 

= Angular coordinate of the wake centerline and wake edge, 
respectively 

= Wake transverse' distances normalized by and L^, respectively 
= Blade stagger 
= Turbulent eddy viscosity 
= Kinetic viscosity 
= Effective kinetic viscosity 
= Density 

= Mass averaged flow coefficient (based on blade tip speed) 


Subscript 

o,e 

c 


t 

2D 


= Values outside the wake 
= Values at the wake centerline 
= Wake averaged values 
= Tip 

= Two dimensional cascade 
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APPENDIX II 
Numerical Analysis 


A. 1 Basie Equations 

Some of the notations and technique used are illustrated in Figure 60. 
Other notations are defined in Appendix I. 

Continuity 


12 + 92 + 32. 0 

ax ay 3z 


Momentum equation in rotating coordinate system. 
X direction (Figure 60) 


ORIGINAL PAGE JB 
OP POOR QUALmr 


1 3p r 3 1 2 3 — ; — ~ , 3 f — p-i 

u' + - 5 — u'v' + — u'w'} 

p 3x 3x 3y 9r 


y direction (Equation 60) 


3V . „ 3V . „ 3V „„„ - 1 3p .3 — T— 7 - ,3 ,2^3 — 7 — r, 

U— +V — + W-r 2f2W cosS - {-r- u'v’ + — v’w'} 

3x 3y 3r p 3y 3 x 3y 9 r 


r direction (Figure 60) 


+ + 2 J 2 V cos 6 - 2 OT sin 6 - rfi 

oX oy oir 


2=_1 

p 3r 3x 3y 


. . .L 9 *2, 

V'W' + T— W' } 

3r 


where x is streamline direction and r is radial direction and y is binormal 
direction. The coordinate frame is chosen as Cartesian at each marching step. 


A. 2 Equations for the Numerical Analysis, Including Turbulence Closure 

The parabolic-ellxptic Navier-Stokes equations at each marching step with x 
as the streamline direction are given below. Lower order terms in the 
equations of motion (Section A1 above) are dropped resulting in the following 
equations. 

X direction: 


„ 3U , „ 3U , „ 3U . ^ 1 3p .3 - 7 - 7 - ^ 9 — t-Ti 

“ ^ 37 ® ■ ■ 7 37 ■ + 8? “ ” > 

y direction: 


„ 9v , „ 3V . „ 3V - . „ 1 3p .3 , 2 3 - 7 - 7 -, 

- - 7 3? - '' ^ ’ 
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r direction 




u||+v|^ + w|^ + 2f2v cos3 - 2OT sing “ ^ 

Here a turbulence model for the Reynolds stresses is introduced. 


3U. an. „ 

-ulul = u^,(7r^ + "^r^) “ ■ • 

X 3 ^ ^ 


Then the momentum equation becomes 
X momentum 


,.au . „ au . „ au . „ i ap . . ,a^u . a^v . a^u . a^j , 

U_ + Y ■—+ \^ — + 2fiW sxnB ~+ y + ■• , H y + 

ax 3y 3r p 3x t\ 2 ax3y 2 3xar 


ay 


ar 


y momentum 

„ av . „ av . „ 3 V . „ 1 ap . . a^v 2 ak , a^v , a^w , 

3x ay 3r ^ p 3y t ^^2 3 3y ^^2 3r8y 

r momentum 


TT „ aw ^ „ 3W ^ „ ^2 1 3p . r3 W ■ 2 3k 

” '5^ H ” w ~ ° ■ F 3^ ■ 3 17 

oT 

2 2 
. 3 V ^ 


■ 3y2 


A. 3 Marching to New x Station with Alternating Directions Implicit Method 

The periodic boundary condition is imposed on boundary surfaces where y = 
constant. Therefore, the first step is r-direction implicit and the second step 
is y-direction implicit. 


A. 3.1 r implicit step 

A. 3. 1.1 Equation for new streamwise velocity components at the new station 

TI.,,.,. - 

2Ay 


U ^j>k + v ^jH-l»k ^j-l-k J_lktl ^I’^l 

j’k Ax ^ 


j’k 


j’k 


2Ar 


+ 2QW. ,, sing 
j’k 


_1 

p 3x 


»T. ■■ >1 


u. 


- 2u. , + u,,. 


= _ i ^ H. V [ i-l’k ^ J’k ^ "1-t-l’k ^i,k-l ^’-‘.i,k ‘ -j’k-fl 


(Ay)' 


(Ar)' 
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4AxAy‘-^j+l’k+l ^j+l’k-1 ^j-l’k+1^ 

4AxAr*^^'^j+l’k+l ^j-l’k-1 “ ^j+l’k-1 “ ^j-l’k+1^^ 

This equation can be rearranged as 

Au Bu . , + Cu . „ , = D 

j’kfl j’k j’k-1 

This system of linear equations with boundary condition can be solved easily with 
standard ADI scheme. See Figure 60 for the definition of U, V, W and u, v, w. 


A. 3. 1.2 Equations for first approximation for v,w in the new x station 

With the upper station value of static pressure, v and w are calculated. The 
finite difference form of v momentum equation is as follows: 

^i’k ^i’k . „ ''^i+l’k ^i-l’k , ,, ^i’k+1 ^j’k-1 

U .,1 ;; + V.,, wr — ^ + W.,. yr — 2flW,, cosg 

3 k An 3 k 2Ay 3 k 2Ar 3 k 


. . 1 Vi-k ' "j-i'k ^ 

p 2Ay t 


Vrk - ^^.rk Vrk 
(Ay)^ 


1 ^3+1'k ~ S*k ~^i’k+l ^~^,rk ^i’k+1 

“ ' (Ar)2 


+ 






4AyAy""3+l’k+l ' j-l’k-1 "j+l’k-l "j-l’k+l^ 

This equation can be rearranged as follows 

A'V.,^_lH-B'V.,3^-bC'V.,^^l = D’ 


Similarly, the equation for w radial velocity at new station can be written and 
rearranged as follows : 




= D" 


A. 3. 1.3 Correction for u,v,w 

Static pressure variation is assumed in the strearai'j’ise direction for the 
parabolic nature of governing equations. Generally, the continuity equation is 
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not satisfxed with the assumed static pressure variation. The static pressure 
is re— estimated so that the velocity components satisfy the continuity equation. 

The standard Possion equation is solved for the cross section at given stage. 

The peculiar point in correcting process is that the mass conservation should be 
considered between periodic boundary surfaces. 

The finite difference form of correction part of static pressure is as follows 




p(Ay ) p(Ar) p(Ay)' 


p(Ar)‘ 
- 




Ax Ay Ar 

and P„ = P. ,, + Pv„ 
j’k j’k j’k 

The equations for correction of velocity components are as follows 

"k 11 

J’k p Ay'-Vl’k 


‘'!’k = -i^3'!-k+i-^!’k) 
"f’k “ " 


and 


e p 

^9xM’k „ 1, , 

S ^AyAr 

P 

where m is the flow rate between periodic boundary surfaces. 


A. 3.2 y- implicit step 

The similar equations to those in 3.1 can be easily developed for y-implicit 
step. At this step of calculation, the periodic boundary condition is applied and 
this step is much like correction step, whereas the r-implicit step is much like 
predict step. 
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A. 4 Calciilation of k and e for Subsequent Step 

The effective kinetic eddy viscosity is calculated front the equation 

'’t - T 

Therefore, at each step k and e should be calculated from the transport equations 
of k and e. 

The finite difference forms of these transport equations is straightforward. 

The equations can be rearranged like 

Ak . , . - + Bk , , , + Ck . , , , — D 

j k-1 j ’ k j k+1 


A. 5 Further Correction for Calculation 

The above sections briefly explain the numerical scheme to be used for the 

✓ 

calculation of rotor wake. The turbulence model includes the effects of rotation 
and streamline curvature. The static pressure distribution will be stored and 
calculated three dimensionally. The result of calculation will be compared 
closely with experimental data. The possibility of correction on turbulence model 
as well as on numerical scheme will be studied. 


ORIGINAL PAGE IS 
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Figure 3 : Curvature of Rotor Wake at r/r^ = 0.721, 1 = 10“ from Rotating Triaxial Probe Data. 
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[3] Rotor wake data, i = 10° (C^ = 0.30) 
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Figure 8: X/C = D. D^2, R/RT=D. 721, I=1D, THETA(DEG„ 1 
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Figure 9: X/C = D. D21. R/RT=D, 721, 1=^10. THETA (DEG. 1 
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Figure 16 : X/C = D. D42, R/RT=0. 721, 1=1 D, THETACDEG. 'i 
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Figure 19 : X/C = Q. D42, R/RT=D. 721, I=1D, THETACDEG. D 
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Figure 21; X / C = D. 02 1 , R / R T= D- 7 2 1 . I = 1 0, TH E T A <DEG. :i 
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Figure 23: X/C = n. U^2, R/RT=D„ 72 1, I = 1 D, TH E T A < D E G . D 
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Figure 25: Static pressure distribution across the wake (x/c = 0.021, R/R_ = 0.721, i = 10°) 
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Figure 26: Static pressure distribution across the wake (x/c = 0 

1 = 10°) 


083, R/R^ = 0.721, 


ts\D" 



CITADEL 0 NO 6<*3 - CI 70 SS SECTION • ZO SQUARES TO INCH 


•iaaiaaiaai ° 

-aaiaaafiaaaia 
aaiaaaBaiaia 


laaaaaa 
laiaaaB 
laiaaBB 

J!HS[ 

S SBia8aaiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiilii| 

aBiaBBBBaaBBiBiaaBaBaaaaBiBBBBailiiaaaaaBBBBBal 
. iBiiiiBBBaaBaiiBiaBBaiaBBiaiBiBBaiaBiaBaiaBBBBBai 
aadpBiiiaBBBBBBiaaBBBBBiaiBlaaBBsBaBaaBBaiaaaaBBaal 

iaiiaiaBiBBBiiaBBiBaBBBBaiBiBaBBaBaBi 

iiiaaliiiBBiiiBRiiaBBBBBatBiaBiiaBiai 
iBaaaaiiaaBiBB8BBiBaBBBBaBaiaaaiaaiBi 

:ii[;s3s»s8s:;s:!;8Es::8;s8t8ul388i 


_iBaaa: 

laaaBBi 

"aaaai' 


aaaBBBBaaBaBaBBaaBaaaBBBBBai 
.aaBaaBaBBBBBaBBaaaaBBBBBaBai 
BBBBaaBBiBiaaBaaaaaaBgBgBBar' 


aaaiBBBBBasB 

“ ~ai: 

iiiiiiiiaaSiiiiiiii 


laBaBaBaaBBaBBBBBaaa 
'~^~aaBKBBiiBBBa 8 aaa 
'BaBaaBiiaaaaaaBB 
naBaBaaBBiaafiaa 
iBataMaBalaBIlaa 



"fBflaiBiBBiBaaiaBBi|BBaaiaBBBBB|BaBai|aii|| 8 !ii 
‘‘"‘aiaiaaiBaaBaBaBBBaiBBaaBBaBBBBaiililalailii 

a 3 ^j===_==j_j_| 


iiiiia 

Hill 
s lllllu 


m 


!S! 


tBBaaBBBB 

-:||BlaB - 

[[•iaf 

BBiiaa 

.aaiiaa 

lissii 




llllllll 

asKssibi 

iaiiaaBBi 

BBBiBBBBI 

tsar"" 
laai 

ail 
ill 


M 


n IH;! 


fBi 


I BBVL^ 

aaiai 

■BiBiBBBBB 

iBiiiaaaa 

jallBaaiB 

iBaBaaaaa 

BBaaaaaaa 



BBBBBBj 

iaaaiBi 

Shb::! 

IIHI 


[lip 

i| 

III 


PI 


U 


.... igiaa 
II 89 «§ 


BaaaBMaaaaaaiBBBBBa 


laiBaaaaaaa 
. iBiaaaaaaai 
IBBBBBBBBBai 


laBBBBaaiBa 

BSBBBBBBSUBI 

lisSSaBBsilhl 


laaaa 

aiaa 

BBB~ 

BBS 

;SB8S 

'ifiK* 

BBaa 

iBaas 

iiU 

S afi- 

r- 


jaaalii- 

aaiBBaaaaBaaaa 
iBBaBBaaiBaaga 
aiaaaaaaagiaii 
Biaaaaaaaiiaaa 
laaaaaaaaa 


gBaaaaaaaaaaaa 

iBBBIBBBBBSBr 

sibsIbbbbbbs 

B BBaaaaaaaBaai 

paaaBr““““““ 

El 


laiaaigaaaaaaiiii 

iBaaaaBaiBaalilii 

EBSBBBliSBBBiBII 

IB! 


S 8«aa 

assB 

BBEr 


IS 

wiiiiaiiiiSB 

lisSiSiE 


aaBB - 

BBEBB 8EBBBE 


SI 


|ai 

ilaa 


.laaaaBMaa 
liaaaBaaaaa 
liaaaaigNaB 
..^^-jiaaiaaEgii 

laaBigiBBgaiaiaa 

BBBBEIIBBBESSBBB 

.BaBaaaaaaaiaiiBa 

aaigaaaaaBaaaaaa 

BaliaaaaaaaaaaEi 

-‘-aaaBBiaaaaaaaia 

BBBBBBBSBBBBBiB 

|aa|aiaiaB|aa|B 

laaBaiBaaBiaaBa 


^aaaaia 

fwsBs: 


—..aaaaBi 

igBaiaaai 

Ubbbbii 

-_.,- 3 ~..jaaaBBaaa 

Issspilslisiili 

lillllllll 

a aiaafigBala 

ishilBBSSSS 


iiipii 

SSSlSi 

■BSEEBBi 


i 8 a 

liaaaaaiaaaBBg 
lEaaiaaaagBaaB 


iaaaai 

Inw 


gaaaaaa 

Baaaaaa 

aaaaaii 

Eisiiu 

m 

as 8 &;s 

aaaala 

: 8 :» 

J 8888 

aaaaaa 

:S88S 

a aaaa 

8IIS 

888 r^ 

■i"< 

3881 . 

aaaai 

B8I8[ 

88888 


i; 


aaaaaiagai 

illsllssi 

IBBSEBBSi N»i 

slSSlijli fiis 

SEBfBsiil 9 *B| 


jliji 
lllllf’ 


18 


IS! 


mill 

iiimi 

IPIHHIH 

II 


Jli 

IB! 


llllirSBI 

lili 


»!p|!U 

SlllSISSalSaBliB 

iBSiliiSSiSSSSi! 


B ebbebi 

maaagi 
— ^aaaBaaaii 

laaiaiaaii 

laagBaaiaaBii 

be3beb8Sb;:i 

SiBsar ■ 

iSi 


JBI 

lil 


III! 

mil 


HI 


laaaaapgi 

ipy 

III 


liiiH 

BBIBBB 

SEbbS: 

■BIS 

ia» 

aaiB 


11 


ail 

III 

BE 


S BBBBSBBB 

ebbSEb:: 

aaaaBaaa 


HMIlilll 

gaaiai jlgi iga 


BBaas 

BSBaSBB 

liamai aa 
1 ! 


il 

fPl 


SS! 


SEE 

nSB 

EEi 

8BE3 


-III 


EEEEbIEb SilE 

m 1!1! 


.aaiiSaaBIS 

aaaaaa 

— aaaa 


iaBBIaai 

aBiBBBaaiiaggiaaaai 
aglaagBliBBigiBaBai 

ips |e I pi i!!l 
lf!||P|!E!||j 

aaaBB&|I||a|Bi i; 

Blilllii 


kaaaaaBHai 

iip 

lijlll 

I 


13 


bEE 

1B8 




a 





IE 





II 


lillllllllll 


■ a a a aa aaaaaaaaaaaaaaaaaaaa 

laaaaBaaaBBBaaBaBBBaaiiiBa 

aBBBaiBBBBaBBaaaBBBBBaaiiaaa 


■BaaaaBBBaaai 


iaaaai BaaBaBaaaaaaaaaaRBaaaaaa. — 

^BBBaaBBaaaaBaaaaBaaaaRaBaaBaaiaaaBBaaBaaBBBaa 

aiii aaaaaa aaaaaa a aaaaaa a aaaaaaaaaaaaaaaa a iaaBBaBBaBBBBB|i 

aiSiBiigSMSSSaaRaaaaaBaaBBBBBiaaaBBBBBBBBaaaBaaBaaaBaBgBaBaBBBaaaaaaa 
aaiaiBi 3 gBiBBBaaaaaaaaaaaaaBaaaaaaRBaiaaaaaaaaaaBaaaBBaBaBBBaMaBaaaaaaa-_---g--------------- 

iaaBiaaaEBBBBBaaaiaBBBaBBaBBaBBaBBRaaaaaaaaaaaaBaBaBaBBBiBBBaBBBaBaBBaBaaaBaaaaiaaiBaaBBBaaB — 

EBiE;EfE!EEEEEEEE:BEBE;EEEEEEEEBEE:!E|E;»^^^ 

aaSBBElSaEBEaBBBBBSBBBBaaSSaSBSSBBaBBB”'"""""'”^”***”””***”*”””- 


SBaSaaBBa 

aaaBBBial 

aaiaBBiaa. 




ESsBEsKB;KKSS|KSS:SS:888S:B8»S3|n 

»»»»»!»esbEbb8s:b::8SSSb:eE8Ebss::sb:|es!3Bh:e8» 


iaaaaaaaaaaai 


'SEBBi 


iBaBaaaaagaaaaaaaBBaBiBaaBi 

jaaaaaiBiiBaaBaaiBaaBRaaaiL 

.aaaiiBaBBBBBBBBBBaaBaaaaaaaaaiaBaaiaiaaaBaBaBBaaBRBaBBBi 
’aaagiBaBBaiaaaaaa aaaaa gaaiaaaaaBaaiaaaaaaBaaaaaBiiBaaaii 

!EEbSS88SEB8S:b:3SBB3S:8S8E88S3:eE;B8SSB8S 

SSBBBEBBSEBaBSSSSgS E!EEbE3SBb!8BS:88!s88S3EEE3BE:SBK 


b:;8B:::ebe3I{bbi;»i 
^”“-”aS3SS!s!BSr** 


i aaaaaaaaaBaaaaaaa 
aaBggBiaalBBBiaaa 
BiBgaiiaaiaaBiBaa 

laiiiaiaaaaaaaaBaB 

aSiSaSaaaaaaaaaaaiBiaiiaiBiiiaaBaiiiiBaiaaan 
i 8 aiia 8 aa|aaaaBaaaBaaaa 3 aaBBavaiaBaaaaaaia|aai 
ilaaaaNiiBiaaBBBBigiaiiaiaaairBiaaMaaiaalaaaaii 
8 gSaaiaag|gaaiaBBiaBaaaaaaBaiaaaBaBaaaalaaaaii 


.aaaaaaaaaiaaaaaafaaBaaaBaBBaBaaaaaaai 

--^aaaaaaiaaaaaaii 

aaaaaE«aaBaBaa'' 


...JEEIbEEEEEeIEEeeEEbEbsEEEEsEEe: IEIHEeb 

iEiEBE-EEEEBBBEEBBEEEiEEEEEEEEEBEE 


aaaBBBBBBaaaiBiaBaaa 

iBaaaBBaaBiBaaaaiaiaaaai 


iiaai 


ialBiBiiiiiaaaaaaa 

'■■■pppBBBaaaBaaiaaaaaar 
:BBBiEiaaaaaBaaBiaBaBaBi 

S gaBaBiaaBBaaBaaaBaaBaaaaaBBBaBaaaaaaaai 
EaaiaiaaBaaBBaaaaBaBBaBaBaBBaaaaaaBBar 


iSEI 

111 

_jaai 

iigagaiaaBaaBBaai 

I igaEalBaaaaBBBai 
liaaBiBgiaEBaSgi 
gBiiiiaiBBaaaaii 

laBBBBBB 

iaaaaaBa: 
iiaaaaaB: 
aaiaiiar 

l■ 888 Kg|■■ 8 BKl 

liE|EEsiB»a{ 

EiBSEEBBlEBBi 


aaaaaiaiaaaaaaaai 
aa aaaaa laaaBaiaaB 
aaaBaiaaaaiaaiaBB 


-jBBBaBaaaiaaaiggaaamBaaBaaBBi aBBBMaaai 
BaiBaaaaiaiaBaiiiBBaaaaaaaBBa 

.gi 


aaaaaaaa , 

aaiaaaaaa a«aaaaaa a aaai 

e 


H| 8 < 

iiai 


aaaaaaaaaaaBaBaaaaai 

— aiaBBiir— 

aaaaaEai 

HKKiBSBBEBBi 


aaaaaaaa 

IBBBiBaaaliaaBBBiaaaaBr''"" 
iaBBiiaiaaaaBBiiiiBiii 

BaBBlBBiaBaaiBiBiiiiEai 

laaBaBiiBaiaaaaaaaiailfBaiaiaii 

^ ^aaaaBi8aBaaalBaiaaaaaaBaiiiiaB||ii|j 

iaaaaaBaaaaaaaiBBaBBggBaaaaaaaiiiB|aBBg|aiii 

aiiaBiaBaaaBBigBaaBBBBBBaaaaBiiiaaBBBaaiaiaiBiiiNgii 
Ba 8 BBBaaBaaaBiaaaaBBaBBaaaBaa 8 aa« 3 BiBaaaaaBia 8 aai 8 i 
■■■■■■BBR»i*a«BMBaBBBaaaaaaai 


aai 

aaaaaaa 

jaaBiaaBBBaaBiBBBi 

iifiliilliiiiil 


'BaaBaag 

iBBauBBa^ 

iiifaaaaai 


aaaaaaaaaai 

— iiaagr***-' 

Baaaaaa. 
aaaaBaai 

I aaaaaa ' 
iaaaar 

8a8888SiaaB888aiaiiaaiiiiiiiiiiiii8iiiiiii8iSi8lii8j 

~aaaaBaaaBaaaaaaBaaaaaaaaiaaiaEaaiaaBaaa*aaaaaaa8aai 

BaaBBBaaaaaaaaBRaaBaBaaaaaaaaaaaaaaai — 

— RaaRRBaaaaaaRBBBaaaaaaaBBaavar-***" 


RRgaRaaai 

Rtla““““ 

fir"" 


fiEssfii 

lliilElE 

liillil! 


aiaa ia 

saiB BSi 

iiaa aaa 

BRii BRB 

Bgii aar 

___ BBBE SSL 

aaa aaaa aa« 


laaaBBBHmHaBRHHaaaBaBBaaaBaa 

llBlpfifiBlifiSfilfigii^lEHfiiBBSSfiBaaiiiiiiBiaRMigMiiiaiiaiiaiiigi 
..B8BBB.i8giEIBik||PiJ8oapaapa|aaa..aa.aa.apai8a»^^^ 


pBPPgpppapPRPPPiPPBaHHaaipSai aip 
= ^iBBBBaiaaiBRBB — 



_^_jBll!BjSBL 

I iaaiaayCSJiiaHaaHaaaaa 
iRaEpaEHiBiaREffippiBP 
ippEEEaalBiEliHapaaap 

EEigPBBBBr" — “““ 

iaiEBaaait 
lagar- 


B8B8:8|SBS|SE| 

bbEbSeEbesEbbEl 

laaaaaBagaaaBBagaaaaiaaBaaaBaBRBa 
laaaaaBiaBBaaiagaaiiiaaaRHBBaBRaa 
jaBiaaaiBBaaaBaEaaiBaBBaHaaaRagBB 
laRPPaaaBHPaPaaaPRBBBapaaaaaRaaaBaaEBBaBlBpaaBaRaa — 

ii88P8PPS8PP999P9!9i!PPP8999!i!!!!!*****"'*'*'*'*9''*** 

■■SSS8SSS8hS!SB 

El 

iBBgMaaaaBMaaiBaBRaiaiiipBBBBiaaaaiiBiiiliBaaaiiBai 
liaBPagaaaBaaaBsaaaBaBiBBBBBBRaaaiPBgEaHaiiaHBl 
i 8 Raaa 8 aaRHRaBRB 8 BaaBaigBBBiBiaaaiBRgiaRagiPRB| 

^BiiiPPpBBRPPapppiiBRpiipRBi 

iBiaERBBBBRaiaBaBaaaHaaaaHB'' 
laiBBigBBBgaiiPBBPBHaaaaaRP 

S iiBiiBBBBaaiaiaaaiBPPiaBa 
BBaiBiBBRaiipiaaBBBaaiaga 
e — — aialp 


aaa aaa 


^ BBaaBBaaBBaaRBi 

8B8SfiSS88SS88!8S88:8SE38fi8LBP iM 

"888888888881888888888:8888 f 


n fiiBBiBa 8 aaRaaaaiiBNiiaaaaapL.^„^...,-„ 

22g"8B8aiRS9iPaS999PPPP9P!il!9!P|99! 


BBRI 

aaSi 


..BBBBRaaai 

BBBBBaaaai 

PBBRBBaaaaaai 

'B8a8aSa8SL 

ggiBBBi^8888888S8BiiiiiBiBiBiitigiaBBBw_.— . — — . 

aaaBBaaBaBBaaBBBBaaaBaaBflaaaBaaiigBBaiaaBa88aiBia|B|BBHa|; 

aaaa aaaaaaBaaaBBaBaBBBBBaaaBBaaaaaBaBaBaiaaaiBRaaaBaaaaaiaiaBHaiaaaaaaaaaai 

Wi IHb iE|l8l8B8!iB S88S888S3 sBB8S38S 888:888:888:88:8888^ 

ml lassi&fiissSiefiiis 88B8s::8888:8888s assfiBBSssSsssssssBS! 

mt 8|888!8fi8SBS88888S888i;B8888::88:S8 88888S8:88:8:&38S8S|Riiiii|MiiiiiM 
IfiS: 8S|iS888888888|EsS888E:88| aiiSipiiiB«Ma8aBa3l,BgBgSI388«Spg8*i8fi^ 

~ iiaaaia aaiai BaaaaBBBBiEBRBaaBBBiaaiiaRPia: 


BaaaaBaa aaa 
BaaRagaa aaa 

— •“-'iPB aaa 


aiBBiiiaaiBBBaaaRBBaaaaaaaaaaaaRaaBaaaBiiaaai 

liBiEiiBaiBaaaaaRaaaaBRBapaaaaaRii laaaaSaaEi 

iiBlgaiagaaaaaaa«apaaaaaBiBaaaaR|| ° 

ifiaaaREiiaaaaiiaRBBBagaigiBBBBBHil 
°-“^-'‘'“aaEaiaag|BaaaaBiaa|aaiaaBRai 


EiE 3 iEsESE 83 EsiSL, 

aBBilalaaaaaaagBRBl 

8ig|aP|aaaBaaiEpgBi 
ipiBapiaiiaaagiBiEi 
„..„aBaal3laiiEaaBEaaai 

BMaaPaaaaaaaiiPaaaPagEBaaaaRaagaaaBaBaBi 


BaaaaaaaaaaaaaaaaL. 
aaaaaaaaaaiapaalSai 
“"-‘“BaaaaaaaEBaaaBi 


SIEiEEIEEiUSEE 


SaaaagiBRaaBaaagaaaaaaaRai 

llaaaaigRaBaaaaRaaaaaaiRai 

SiiBaaiiaai8PPaapaaaiiiaa| 

nBBBiigBBiBaaaaiiHBl 

liBaaBBBBEaaiar 


jia!Hinjl|Pl|: 3 ESaE 3 H 8 |asEfj» 
IfiEsSElisiSSsIlEEEEEESiaEsEEiiEsHai 


SE;;iE:sK;»isu 

— ==-paB|a*aB|-“ 

"ifii 


aBaaBaaag 

-aaaaiiii 

aaaaaEia 


sliSEKSSai 
ifii4B38S|:8EEL.. 

JiaiaaRaaaEiBaaaaai 



ar' 


mi\\\ Iplilpp! 

laaaglaaaiBE iiBaaaaBiai 

... ... 

*'li 


aBBaaBiiiaa Biaaai . 
aiiiaiiaaBi aiaai bl.. 
aalBBaiaBBa iaaai iaai 

, EIESEEESEEE SEISE EEEi 

Hilililimiimiilliilllllliilif 


,ise EE! 

EtiEi aia 


„aaaaa_ 

laaaaiRapal 
>aaaaa«aal 
EaRaaaaaal 
BaaaaaaaR 
aaBBaatai 
aaaBaaBBl. 

iiaaaaaaal 


118881811 


jaaRaa 

aaaBaaiaE. 
naBBaBaail 

B aaaaRail 
aaaaRaal 
laaagaRaal 
— laaRi — ■ 


aaBaBBaaal 

JEESEEiESl 


aaaagaaai. 

aaBBiBHaal 

aaagBBBaal 

SEe3s8|8sI 
eeeeeSiesI 

ESEEESar 

lagaaaaa 
laSaaaaa 

iiiiiiiiiii 


sIIIIeIIII 

BaaaaaSSfil 

■agaBaiSfi^ 

SSKSESEi 


EES 

PBR 


;i 


. aaa iaaaag 

*BBBag 

aiaai 


lliiSbshllSsflhsil 


>88S8S888888888888S8B 

“'"“^“aapaiaaiaiigpp. 

aapaEaBa8j|Sii|| 


!g l 

I 
I 

BBaaaBaaBBL^ 

E SEi!E!E!EI!EE|EEE|sE 

BEIlilEiaiiEnilEHl 


Figure 28: Static pressure distribution across the wake (x/c = 0,021, R/r^ = 0.860, i = 10°)^ V0‘ 
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FIG. 54 FOURIER COEFFICIENT VERSES Z/C 



FOURIER COEPFICIENT 

- 0.2 0.0 



0.0 


0.2 


0.4 


T 


0.6 0.8 

NORMALISED AXIAL DISTAl^CE Z/C 


1.0 


FIG. 55 PLOT OF FOURIER COEFFICIENT VERSES 2/C 
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AXXAL MEASURED UD/UDC, FOURIER UD/UDC 

Kgure 57; Goraparison o£ Fourier Curve and Data 
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